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We characterize phases of the compass ladder model by using degenerate perturbation theory, symmetry 
fractionalization, and numerical techniques. Through degenerate perturbation theory we obtain an effective 
Hamiltonian for each phase of the model, and show that a cluster model and the Ising model encapsulate the 
nature of all phases. In particular, the cluster phase has a symmetry-protected topological order, protected by a 
specific Z 2 X Z 2 symmetry, and the Ising phase has a Z 2 -symmetry-breaking order characterized by a local order 
parameter expressed by the magnetization exponent 0.12 ± 0.01. The symmetry-protected topological phases 
inherit all properties of the cluster phases, although we show analytically and numerically that they belong to 
different classes. In addition, we study the one-dimensional quantum compass model, which naturally emerges 
from the compass ladder, and show that a partial symmetry breaking occurs upon quantum phase transition. We 
numerically demonstrate that a local order parameter accurately determines the quantum critical point and its 
corresponding universality class. 
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I. INTRODUCTION 

A comprehensive understanding of the phases of mat¬ 
ter and also type of transition between them have long 
been a principal problems in condensed matter physics. It 
was believed that Landau-Ginzburg theory can help provide 
such understanding.^ This theory is based on ‘spontaneous 
symmetry-breaking phenomenon’ associated with a nonzero 
‘local order parameter.’ In this theory, all phases of matter 
are identified by some broken symmetries or, equivalently, 
by their corresponding local order parameters. However, the 
emergence of the so-called ‘topological phases,’ which has no 
evidence of symmetry breaking, defies this theory.Topo¬ 
logical phases manifest exotic properties such as robustness 
against local perturbations^’®, nontrivial anyonic statistics^’*, 
and exhibiting long-range entanglement®, which make them 
interesting theoretically and experimentally. 

In the past two decades, vast efforts have been devoted to 
providing ‘an alternative framework’ for characterizing ex¬ 
otic phases of matter. Recently, inspired by ideas from quan¬ 
tum information theory (especially distribution of entangle¬ 
ment), “symmetry fractionalization” has been proposed as 
a technique for full classification of the phases of (quasi) 
one-dimensional (ID) gapped quantum systems has been 
proposed.*®"'* This classification, based on structure of entan¬ 
glement, places the phases into three classes: (i) symmetry- 
protected topological (SPT) phases, which have short-range 
entanglement, (ii) topologically-trivial phases, which can be 
mapped to fully-product states (with zero entanglement), 
and (iii) symmetry-breaking phases (with degenerate ground 
states). SPT phases, unlike topologically-trivial phases, can¬ 
not be mapped to a fully-product state as long as some specific 
symmetries are preserved; that is, they are robust against any 
perturbations which respect these symmetries. 

In symmetry fractionalization, one needs to determine those 
symmetries which protect a phase, from which a set of unique 
labels are obtained to distinguish the phases that are sepa¬ 
rated by a quantum phase transition—see Sec. IV C. Obtain¬ 


ing phase labels, however, is a challenging task, which gener¬ 
ally requires the prior knowledge of symmetries of the model 
and also an exact infinite matrix product state (iMPS) repre¬ 
sentation of its ground state. Having determined the symme¬ 
tries and the iMPS representation of ground state, e.g., by us¬ 
ing the infinite time evolving block decimation (iTEBD) or 
infinite-size density matrix renormalization group (iDMRG) 
methods,one can employ the techniques proposed in 
Refs. 16 and 17 to determine phase labels. 

There exist numerous (exotic) models which have been 
proven to exhibit topological order, but yet a simple and ex¬ 
perimentally realizable model featuring topological phases is 
of great interest.**"*® In this respect, the Kitaev honeycomb 
model has been a prominent candidate.**"*® The Hamiltonian 
of this model contains two-body interactions (hence relatively 
easier to realize experimentally), and has a rich phase diagram 
that exhibits different classes of topological phases and non- 
Abelian anyons. In addition, the Kitaev honeycomb model 
on an arbitrary-row brick-wall lattice (another representation 
of the honeycomb lattice) has also been recently studied.** 
The associated quantum phase transition between the ‘exotic 
phases’ of these models are believed to be of topological type, 
without any (spontaneous) symmetry braking. Nevertheless, 
the characterization of these phases had remained largely un¬ 
known; this is indeed our very goal here to bridge this gap. 
The model on one- and two-row brick-wall lattices takes a 
simple form referred to as the “ID compass”** and the “com¬ 
pass ladder” models, respectively. Characterization of the 
corresponding phases is of special importance because these 
phases (with a proper modification) also appear in the phase 
diagram of the Kitaev honeycomb model on arbitrary-row 
brick-wall lattices. In addition, since ladder systems can be 
created and manipulated by highly-controlled quantum simu¬ 
lators, they play an important role in experimental realization 
of ‘Majorana fermions’—and whence topological quantum 
computation.*®"** A promising platform based on the ‘inho¬ 
mogeneous Kitaev ladder model’ has been recently proposed, 
which can read out Majorana fermion qubit states and also 
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perform non-Abelian braiding. 

Our main objective in this paper is to identify the type of 
quantum phase transitions and different topological phases of 
the compass ladder and ID compass models. The compass 
ladder includes three phases denoted by 21, 25, and £—see 
Fig. 1. We employ degenerate perturbation theory,^^ to as¬ 
sign an effective Hamiltonian for each phase, which yields: 
(i) (two different) cluster model(s)^^~^'—written in different 
basis—for the 21 and £ phases, and (ii) the Ising model for the 
*B phase. Based on this analysis, it is shown that the 21 and 
£ phases belong to the cluster phase, which is a well-known 
SPT phase protected by the Z 2 x Z 2 symmetry. Despite sim¬ 
ilarity of the 21 and £ phases, we show that they belong to 
different classes of SPT phase; the 21 phase is protected by 
the complex-conjugate symmetry, whereas the £ phase is not. 
This observation is also numerically verified by the iTEBD 
method and the symmetry fractionalization technique. 

The 25 phase appears to be of topologically-trivial Z 2 - 
symmetry-breaking type, characterized by a Landau-type lo¬ 
cal order parameter. This implies a spontaneous symmetry 
breaking upon quantum phase transitions, and thus, the phase 
diagram of the compass ladder can be classified by the asso¬ 
ciated local order parameter. We demonstrate this result by 
the iTEBD method after determining the local order param¬ 
eter and symmetry-breaking group—see Eig. 4. In addition, 
the local order parameter correctly specifies the universality 
class of the quantum phase transitions as of the Ising class 
(with the magnetization exponent f3 = 1/8). We remark that 
our conclusion differs with Ref. 27, where the classification 
of the phase diagram is based on nonlocal string order param¬ 
eters (whereby believed that there were no explicit change of 
symmetry upon quantum phase transitions). 

Additionally, we study the ID compass model, which natu¬ 
rally appears by turning off one of the coupling parameters of 
the compass ladder. Upon quantum phase transition a spe¬ 
cific Z 2 symmetry is broken and another one is preserved, 
thus a partial spontaneous symmetry breaking occurs —i.e. a 
quantum phase transition between two phases, where in each 
phase, part of symmetry group has been broken. Based on 
this fact, one can construct a local order parameter to cap¬ 
ture quantum the phase transitions and relevant physics of the 
model. Interestingly, as examined by the iTEBD method, this 
local order parameter is shown to give the accurate values of 
both critical point and magnetization exponent (/3 = 1/8). 

This paper is organized as follows. In Sec. II the models 
and their phase diagram are reviewed. In Sec. Ill the effec¬ 
tive Hamiltonian of the compass ladder is obtained. Broken 
symmetry of the *8 phase and its corresponding local order 
parameter are derived in Sec. IV, and numerically examined. 
The implementation of the symmetry fractionalization tech¬ 
nique to obtain the labels of the SPT phases are presented in 
Sec. IV C, and the topological properties of the SPT phases are 
discussed next in Sec. IV D. We discuss the phase characteri¬ 
zation of the compass model in Sec. V. The paper is concluded 
in Sec. VI with a summary of our results. 
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FIG. 1. (Color online). Graphical representations and phase dia¬ 
grams of the compass ladder and ID compass models. Black circles 
and colored links represent spin-1/2 particles and different types of 
interactions, respectively, (a) The compass ladder model, (b) The ID 
compass model, obtained by switching off the red-link interactions 
of the compass ladder in (a), (c) The phase diagram of the compass 
ladder model. Here the two paths Ui and 32 are introduced for our 
numerical analysis, (d) The phase diagram of ID compass model. 


II. COMPASS LADDER MODEL 

The compass ladder model (also referred to as the XYZ 
compass model"^^) is defined on a ladder geometry as in Eig. 1 - 
(a), where the black circles denote spin-1/2 particles, and the 
colored links (blue, red, and violet) represent different types 
of interaction denoted, respectively, by ‘b,’ ‘v,’ and ‘r.’ The 
Hamiltonian is given by 

HKi. = -Jh X! 

b links r links v links 

( 1 ) 

where (t“ (for a G {x,y,z}) represents the a Pauli matrix, 
and Ja (for a G {r, v, b}) is the coupling constant. Without 
loss of generality, the coupling constants are assumed to be 
positive; Ja ^ 0. In Ref. 27 the phase diagram of the model 
has been obtained as in Eig. l-(c) through the Jordan-Wigner 
transformation technique. This diagram contains three gapped 
phases labelled by 21,58, and £. The 21 (25) phase is separated 
from the 25 (£) phase by the gapless line Jr/Jv = Jh/Jv + 1 
{Jr/ Jv = Jh/Jv — !)■ The quantum phase transition between 
these phases is of the second-order type (because of the diver¬ 
gence in the second derivative of the ground-state energy), and 
was believed to be topological (characterized by string order 
parameters). 

The compass ladder model reduces to the ID compass 
model when one of the coupling constants vanishes. Eor the 
case of Jr = 0, as shown in Fig. l-(b), the Hamiltonian re¬ 
duces to 

iTcompass = - Jh ^ cr/aj - ^ ct/ctJ. (2) 

b links V links 

The phase diagram of the ID compass model is already known 
as in Fig. l-(d), and contains two gapped phases with exten¬ 
sive degeneracy separated at the critical point Jv /Jh = 1- The 
gapped phases are called the blue- and violet-compass phases 
for Jv/Jb < 1 and Jv/Jb > 1, respectively. The quantum 











3 


phase transition is topological and of the second-order type. 
Similar to the compass ladder, the nature of the topological 
quantum phase transition in the ID compass model has been 
shown through nonlocal string order parameters. 

To identify the nature of the 21, 05, and £ phases and their 
corresponding quantum phase transitions, we derive the effec¬ 
tive Hamiltonian"^^ of the compass ladder associated with each 
phase, which capture main physical properties of each phase. 


III. DEGENERATE PERTURBATION THEORY AND 
EFFECTIVE HAMILTONIANS 


Degenerate perturbation theory is based on splitting the 
original Hamiltonian H into two parts: Hq and V. The Hq 
term represents the unperturbed Hamiltonian, whose energy 
spectrum is fully known, and in general could be degener¬ 
ate. The V term plays the role of perturbation, whose operator 
norm is relatively smaller than the spectral gap Aq of Hq, i.e., 
||U|j ^ Aq. In the case of the degenerate perturbation for¬ 
malism, for a specific energy level, the correction at the mth 
order of perturbation is given by an ‘effective Hamiltonian.’ 
For a quantum phase transition, the effective Hamiltonian for 
ground-state energy is required, which is denoted by H^^'^ 

The starting point to obtain is to define the projection 
operator into the ‘unperturbed degenerate ground space’ (set 
of all ground states of Hq), 

‘p= E l^o)('I^ol, (3) 


where Eq is the ground-state energy of Hq. Having deter¬ 
mined iP, the effective Hamiltonian H^'^ can be determined. 
The first-order effective Hamiltonian H^ has the following 
form: 


= ^PU^P. (4) 


The form of higher orders of the effective Hamiltonian be¬ 
comes gradually more complex,; e.g., the second- and third- 
order effective Hamiltonians are given by 

^ rpYGV^, (5) 

H^J^^ = (PVGVGVE - e'^q'^'EVGGV'E, (6) 


where 


G = 


1 

Eq-Hq 


(l-tP) 


(7) 


is the Green’s function, and e'^'’ denotes the ground-state en¬ 
ergy of 


A. Effective Hamiltonian associated with the 23 phase 

To obtain the effective Hamiltonian for the iB phase, Hq 
and V are set as follows: 

Hq =Jv E 

V links 

V=Jq E + E 

b links r links 

where ,Jq J„ (note that positivity of Jv and nonzero 
values of Jj. and Jq guarantee that the ground state of H = 
Hq -f U is within the 25 phase (see Fig. l-(c)). 

The projection operator tPv, which comes from the unper¬ 
turbed degenerate ground space (set of all highly-degenerate 
ground states of Hq), is defined as follows: 

n =itt)(tti+iu)(ui, (8) 

V links 

where | f) and | j,) are the eigenstates of —index v denotes 
violet. However, for simplicity it is more convenient to write 
HeS in a new basis by rewriting fpQ as 

2’“ = |t)(ttl + II)(U|, (9) 

where jt) = | tt) and |j,) = | j,4,) are the ‘logical qubits’ in 
the tr^-basis. The energy and degeneracy of an unperturbed 
ground state are, respectively, equal to Eq = —N^J^ and 
2^'', where is number of the violet links. The first ex¬ 
citation of Hq has energy Ei = — (iVv — 2) Jv with degener¬ 
acy which is obtained by flipping one of the spins. 

Flipping two spins on different violet links gives rise to higher 
exited states that has the energy E 2 = —{Ny — 4) Jv with de¬ 
generacy 2{Ny — l)Av2^''“^. 

The first-order effective Hamiltonian {fPVfP) is zero be¬ 
cause V excites the unperturbed ground space into the second- 
excited subspace, which obviously has no overlap with the 
unperturbed ground-state subspace; whence H^ — 0. How¬ 
ever, the second-order effective Hamiltonian is nonzero, re¬ 
sulting in both nontrivial and trivial terms, (nontrivial terms 
break the highly-degenerate ground-state subspace, while the 
trivial terms do not). In the expression fPVGVfP, one of 
the possibilities (among many) is to choose the first and sec¬ 
ond U on a specific link. The first V excites unperturbed 
ground space, bringing it to the second excited space. The 
effect of the Green’s function on the second excited space is 
G = —l/(4Jv), and the second V takes the second excited 
space back to the unperturbed ground space. Thus, one can 
show that such interactions yield trivial contributions to the 
second-order effective Hamiltonian 



where / = |t) (t| + |i) (i| and i runs over the violet links, see 
Fig. 2-(a) [left]. 

The other nonzero contributions to the second-order effec- 

( 2 ) 

tive Hamiltonian H^g are from those interactions that act on 
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FIG. 2. (Color online). Schematic representation of the effective 
Hamiltonians of the 21, 58, and £ phases. The vertical violet ellipses 
show that the state of the two spins (denoted by the black circles 
placed within the ellipses) is either | tt) or | 44-) ■ The horizontal 
blue and red ellipses also represent the states {| /'/')■, \ i/i/)} 
and {| W), I W)}, respectively (see the main text). Colored 
links and circles denote different types of intersections and logical 
qubits. Index i labels the ellipses [left], or equivalently logical quhits 
[right], (a), (b), and (c) [left]: The main contributions to the effective 
Hamiltonians of the 58, £, and 521 phases, respectively, and [right] 
their corresponding effective Hamiltonians. 


nearest neighbor violet links, as sketched in Fig. 2-(a) [left]. In 
this case, the hrst V —the blue link in Fig. 2-(a) [left]—excites 
the unperturbed ground space, resulting in the second excited 
space. The action of the Green’s function on this second ex¬ 
cited space is given by G = — l/(4Jv). The second V —the 
red link in Fig. 2-(a) [left]—takes the second excited space 
back into the unperturbed ground space. It is straightforward 
to show that ^’VGV (P for Fig. 2-(a) [left] is proportional to 

where is the y Pauli matrix 
in the logical basis {|t), |4)}, i.e., 

a« = -*|t)d|+z|4)(t|. (11) 


We note that the terms acting on the next-nearest-neighbor 

violet links (or farther neighbors) play no role in the second- 

order effective Hamiltonian. 

( 2 ) 

In summary, is given by, as shown in Fig. 2-(a) [right], 




'Jv. 


,Tb Jr 




= -2 


4J, 
Jb Jr 


4Jv 


4Jv 


erf crf+i -F const. 


'i+1 


( 12 ) 


The factor 2 indicates that there are two possibilities for 
choosing the first and second V in the expression [PVGV 5P. 


B. Effective Hamiltonian associated with the £ phase 


where Jv, Jb <C Jr- This sort of definition of Hq, V, and the 
coupling constants is to guarantee that the ground state (of H) 
is placed within the £ phase. Similar to Sec. Ill A, the goal is 
to obtain the leading-order nontrivial effective Hamiltonian. 

The projection operator into the highly-degenerate ground 
state of iJo is given as follows: 

5Pb= n 5Pb°^idd)(ddi+idd)(ddi, 

b links 

where | /') and | \) are the eigenstates of . Rewriting 
in a new basis makes the form of the effective Hamiltonian 
simpler as 

= l7)(dd I + i7)(dd I, (15) 

where \/^) = | /'/') and |^,/) = | v'v') th® logical 
qubits in the -basis. The energy and number of degeneracy 
of the unperturbed ground space are the same as Sec. Ill A, 
i.e., Eq = —NbJb and 2^*=, where is number of the blue 
links. The hrst (second) unperturbed excited space is obtained 
by Hipping one (two) spin(s) on a specihe (two different) blue 
link(s), which give Ei = —{N^ — 2) Jb ( £'2 = — — 4) Jb) 
with degeneracy 2Wb2^‘’“^ (2(£b — l)A^b2^‘’“^). 

Similar to Sec. Ill A, the Hrst-order effective Hamiltonian 
is zero: = Jbl^Jb = 0. The second-order effective 

Hamiltonian results in trivial terms: the only possibility to 
have nonzero terms for 1P14 GV IP is to choose the Hrst and the 
second y on a speciHc link. It yields 

i i 

where I = \/') {/'\ -F |\) (\|, and i runs over logical qubits, 
as shown in Fig. 2-(b). The third-order effective Hamiltonian, 
E[^ = ¥hVGVGV¥h — El^^(PhVGG(Ph, leads to a nontriv¬ 
ial term. The second term of El^ vanishes because E^'^ = 0. 
The closed form of the Hrst term (5PbT^ GV GV 5Pb) is obtained 
by such choices as depicted in Fig. 2-(b) [left]. Suppose the 
Hrst and the second V are the violet-link interactions, and the 
third V is the red-link one. The Hrst V excites the unper¬ 
turbed ground space to the second excited space. The effect 
of the Green’s function G on the the second excited space is 
G = — l/(4Jb). The second V just transforms the second 
excited state to itself; that is, the second V only rotates the 
states within the second excited space. Thus, when the next G 
is applied, G = — l/(4Jb). The third V takes the second ex¬ 
cited state back into the unperturbed ground space. It can be 
shown that the expression ^P^VGVGV^Ph, in Fig- 2-(b) [left], 
is proportional to • • • 0 / 0 ct® 0 ct* 0 ^ j ^ ^ where 

and are the x and z Pauli matrices in the logical basis 
{\/^), li/)}. Here and are given by 


Here Pfg and V are deHned as follows: 

£o = Jb ^ 

b links 

V = j,Y 77 Y 77- 

r links v links 


-\AAV y = \A{A + \A{A- (17) 

Other selections of ‘BaVGVGV^Ph —except those in Fig. 2- 
(b) [left]—make no contribution to , whence 

J P 


777 = -2 


(4Jb) 


2 '^i-E 


'i+l- 


(14) 


(18) 
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The factor 2 is again due to different choices of V —there are 
6 different configurations, similar to that of Fig. 2-(b) [left], 
whose factors cancel out each other as —2 + 2 — 2 = —2. 
Equation (18) is the cluster Hamiltonian, which belongs to 
the class of stabilizer Hamiltonians. The ground state of the 
cluster Hamiltonian has a unique (for periodic boundary con¬ 
dition) and exact MPS form, and is of the Z 2 x Z 2 SPT type."^*^ 


C. Effective Hamiltonian associated with the 21 phase 


The effective Hamiltonian of the 21 phase can be obtained 
by replacing Jj. —>■ Jb and —>■ a^(cr^) in the results of 

Sec. Ill B, which yields 




(4J.) 




(Ta a. 


i+l ’ 


(19) 







ul 


FIG. 3. (Color online). Diagrammatic representation of IMPS and 
some related quantities, (a) Graphical representation of Eq. (21) 
with the matrices A), (b) Conditions for canonical IMPS, 

(c) Expectation value of an on-site operator 0 in the canonical IMPS 
form, (d) Eigenvalue equation of the (^-transfer matrix with its max¬ 
imum eigenvalue (A^) and the corresponding right eigenstate (DJ). 
If A^ = 1, the IMPS is symmetric under tj. 


where and ct* are the y and z Pauli matrices in the logical 
basis {|\),|\)}- In this basis, 

l\) = l\\), l\) = l\\), (20) 

where | \) and | \) are the eigenvectors of Equa¬ 
tion (19) is the cluster Hamiltonian written in a different ba¬ 
sis; it can be obtained from Eq. (18) by 7r/2-rotation about 
the z-axis. Since this operation is unitary, the ground state 
of the Hamiltonian (19) inherits the properties of the cluster 
phase such as having unique exact MPS form and being of the 
Z 2 X Z 2 SPT type. 


IV. CHARACTERIZATION OE DIEFERENT PHASES 

A. Infinite matrix product state (IMPS) method 

Ground state of (quasi) ID gapped quantum systems re¬ 
spects ‘area law,’ in the sense that bipartite entanglement of an 
arbitrary subsystem depends on its boundary rather than bulk. 
Based on this fact, it has been proven that (quasi) ID gapped 
quantum phases can be faithfully represented by iMPSs.^^^ 
The IMPS representation of a state |'k) (ground state of a ID 
gapped system) is based on assigning to each site a set of ma¬ 
trices as 

l^)= I] 

■■■mi.mi+i"- 

( 21 ) 

where A is a D x D diagonal matrix, and E^^'^s are some 
D X D matrices assigned to site i [Fig. 3-(a)]. The matrices 
(rC™*), A) are usually determined by the iTEBD or iDMRG 
methods, where the accuracy of the scheme is controlled by 
the parameter D. Having determined the matrices A), 

one can always use a ‘canonical transformation’ and rewrite 
the iMPS representation in a more suitable canonical form: 
(E^™-), A) -w A)4®. In the canonical IMPS form, as 


shown in Fig. 3-(b), new matrices A) satisfy the fol¬ 

lowing conditions: 

^(f(rn.)^)(p(m,)^)t ^ (22) 

rrii 

^(Af(™-))t(Af(™*)) = 1. (23) 

rrii 

where A is a positive diagonal matrix related to the density 
matrix of a half of the system through g = A^. In this form, 
the expectation value of a local order parameter (defined on a 
given site) is given by 

(vk|d|vi/)= 5] (A(„)rj:-)^)A(^))d(™.),(^,) 

mim'-oc(3 

X (A(a)f!:;i)A(^))*, (24) 

as depicted in Fig. 3-(c). In addition, the on-site symmetry 

groups can be evaluated in a straightforward manner in the 

canonical iMPS representation of the ground state I'F). The 
on-site symmetry Q = 9i respected by IaF) if in the 
following relation Xg becomes 1: 

= (25) 

where Tg is the ^-transfer matrix 

= E (ri7^^/5)(5(2n,).«))(fri^V)*, 

rriim'- 

(26) 

shown in Fig. 3-(d), and A^ is its maximum eigenvalue. Fur¬ 
thermore, if the symmetry Q is respected (that is, A^ = 1), 
the following relation should be satisfied: 

E5(m.).(™')r('"-) = (27) 
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where e*®® is a phase, and Ug is a unitary matrix (which plays 
an important role in the classification of SPT phases)—see 
Sec. IV C. It is straightforward to show that the right eigen¬ 
state of the ^-transfer matrix Tg (corresponding to the eigen¬ 
value A^) is Wg (see Fig. 3-(d)). 

B. Local order parameter 

The nature of the 55 phase is revealed by the Ising Hamil¬ 
tonian (12). This Hamiltonian has two fully-product de¬ 
generate ground states, implying that the 55 phase is of the 
topologically-trivial Z 2 -symmetry-breaking type. The Z 2 - 
symmetry-broken group and the corresponding local order pa¬ 
rameter (O), in the logical basis, are given by 

Z2 = {X,1}, 0 = Y,ayN^, (28) 

where X = /j. By employing the projection 

operator tPv, these two quantities can be recast in the original 
basis as follows: 

o = ^ c» = E )/2^v, 

2=1 V links 

X = \lat^X= n afaj, 

2 V links 

l=\{h^I= J] IXg. (29) 

2 V links 

The broken symmetry group Z 2 = {X, 1} and the local order 
parameter O uniquely characterize the S phase in the sense 
that in this phase the symmetry X is not preserved, and the 
local order parameter O is nonzero. 

On the other hand, the 21 and £ phases represent nonde¬ 
generate ground states, which both respect all symmetries, in¬ 
cluding X. This yields that O is always zero within the 21 and 
£ phases, 

{0) = {og) = -{go) = -{0), g= W 

V links 

^(O) = -(O)=0, 

where the operator g is one of the symmetries of the model. 
As a result, the phase diagram of the compass ladder can be 
classified by the local order parameter O. 

We have numerically plotted the local order parameter O 
through the paths {3i,T2} in Fig- 4-(a). The plot indicates 
that whenever the IB phase appears (in the range of 1 < 
Jh! Jv < 3 and 0 < Jb/Jv < 2, respectively, for paths 3i 
and 32 ) the local order parameter O becomes nonzero. In ad¬ 
dition, O decays when it approaches the boundaries of the 55 
phase—the points Jh/J^ £ {1)3} and Jb/Jv £ {0)2}, re¬ 
spectively, for the paths 3i and 32- As plotted in Fig. 4-(b), O 
vanishes as 

0~|Jb/Jv-3|^ (30) 

0~|Jb/Jv-2|^ (31) 



FIG. 4. (Color online). Local order parameter, its scaling close to the 
critical points and the broken symmetry along the paths {3i, ^ 2 }. (a) 
The local order parameter 0 becomes nonzero within the 23 phase, 
which indicates spontaneous symmetry breaking, (b) Log-log plot 
of the local order parameter 0 versus J^/ X in the vicinity of the 
critical point, which gives the magnetization exponent /3 = 0.12 ± 
0.01. (c) Ax shows qualitative behavior of the symmetry Af; when 
Ax < 1 it implies that X is broken. In the 23 phase, as expected, the 
symmetry X is spontaneously broken. 


in the vicinity of the boundary points 3 and 2 for the paths 
{3i,32}, respectively, where /3 = 0.12 ± 0.01. This implies 
that the exponent /3 is 1 /8, and the quantum phase transition is 
of the second-order type. The same results have been obtained 
by using nonlocal string order parameters in Ref. 27. 

The behavior of the symmetry X can be explicitly investi¬ 
gated by calculating the maximum eigenvalue of AC-transfer 
matrix (i.e.. Ax), as plotted along the paths {Ti, II 2 } in Fig. 4- 
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(c). Again, whenever the 58 phase appears, Ax becomes < 1, 
implying that the symmetry has been broken. This observa¬ 
tion agrees with the effective Hamiltonian (12). 


C. Symmetry fractionalization 


where is a phase. Taking complex conjugate of Eq. (34) 
and iterating this equation twice gives 

r* = ulsu^ ^ r =ul^*r*u*^ ^ r = 


The technique of symmetry fractionalization provides a 
method to uniquely distinguish different SPT phases. This 
technique for ID gapped systems is complete, and provides a 
set of unique labels assigned to each SPT phase. These labels 
are obtained by transformation of the iMPS representation un¬ 
der the symmetries of system. To clarify how these symme¬ 
tries result in unique labels, we shall discuss two examples; 
Z 2 X Z 2 and symmetries. 

Assume that the on-site symmetries Q = Y\i9i ^ — 
hi commute; gihj = hjgi, and = hf = t (for all i 
and j). These symmetries are isomorphic to the Z 2 symmetry 
group in the form of /} and {Q^ /}. One can combine 
these Z 2 symmetry groups and form a Z 2 x Z 2 group with 
elements Q9-[, 7}. If Z 2 x Z 2 is respected by the 

iMPS, the maximum eigenvalue of Q- and fTf-transfer ma¬ 
trices should be equal to one {\g = = 1) and Eq. (27) 

should be satisfied for the elements of the symmetry group. 
Equation (27) yields 

ghf = UlulfUhUg, hgf = ululfUgUh, 

^UgUh = e^^^’^UhUg, (32) 

ggf = f = UlulfUgUg, hhf = f = ululfUhUh, 
^UgUg=e^^-, UhUh = e^^\ (33) 


where the phase factor is used to classify SPT phases 
(note that for simplicity the summations and phase (e*®) have 
been ignored). By Eqs. (32) and (33), e*^*''* can only be 
±1. This allows two different orders: the SPT phase with 
gzOgh _ jjjg trivial phase with = -fl. Through¬ 

out the SPT (trivial) phase, we have e**^9'* = —1(-|-1); the sign 
changes only upon a quantum phase transition. The minus 
sign also reveals that the SPT phase is protected by Z 2 x Z 2 
symmetry; i.e., any perturbation which respects the symme¬ 
try cannot destroy the SPT phase. The two signs also repre¬ 
sent two inequivalent projective representations of the Z 2 x Z 2 
symmetry—see also Refs. 37 and 47. 

Based on this observation, the topological order parameter 
C^ 2 xZ 2 introduced as follows: 


{l/D)Tr[UgUhUlul] 


|A^| < lor|A^| < 1 
= 1 


This order parameter only takes values {0,1, —1}, from which 
the phase can be characterized. Specifically, the values 
0, 1, and —1, respectively, denote the symmetry-breaking, 
topologically-trivial, and SPT phases—corresponding to the 
Z 2 X Z 2 symmetry. 

If the iMPS is symmetric under the complex conjugate sym¬ 
metry = 1 and Eq. (27) becomes 


j* 


Sacf/t 


•)[/, 




(34) 


(for simplicity index {rrii) and the arbitrary phase have 
been ignored). Since is unitary, the phase becomes 
±1. Each of these signs denote a separate order. Specifically, 
gj^x = — 1 indicates an SPT phase protected by whereas 
gj^x = 1 indicates a topologically-trivial phase. Similar to 
OZ 2 XZ 2 ’ 0 * 1 ® define a topological order parameter (Oj^;) 
that detects topological properties of the SPT phase protected 
by K, 


n - f ® ’ I^A;I < 1 

- I (l/D)Tr[Ug^U*^} ; [Aj^l = 1 ' 

D. Topological order parameter 

The 21 and € phases have SPT orders, as we showed by our 
degenerate perturbation analysis. In this section, we investi¬ 
gate the topological aspects of theses phases, namely: (i) there 
is a specific Z 2 x Z 2 symmetry which protects both phases, 
and (ii) the complex-conjugate symmetry protects only the 21 
phase, which indicates that the 21 and £ phases belong to dif¬ 
ferent classes of SPT phases. 

The 21 phase is characterized by the cluster Hamiltonian 
(19), and its ground state belongs to an SPT phase protected 
by the following Z 2 x Z 2 symmetry group^®’"^*’ (written in the 
logical basis): 


Z2xz2 = {g,^,g^,i}, ( 35 ) 

g^l[^2^, (36) 

2i 2i+l 

Rewriting this symmetry group in the original basis results in 
g= W cTig, W (37) 

V links V links 

and Z 2 X Z 2 = {^, 9 -C, g^, 7}. Thus, the associated topo¬ 
logical order parameter C^ 2 xZ 2 should take the value — 1 
(which signals the existence of SPT phase) within the whole 
region of the 21 phase. 

It is straightforward to see that the Z 2 x Z 2 symmetry group 
of the £ phase has the exact form of Eq. (37). Thus, one con¬ 
cludes that 0 ^ 2 x 22 should be equal to —1 for both 21 and £ 
phases, indicating the SPT phase protected by Z 2 x Z 2 ; and 0 
for the 58 phase, implying the symmetry-breaking phase. 

The topological order parameter 0^2 xZ 2 has been plotted in 
Pig. 5 for the paths {3i,32}- This plot confirms that 0^2x22 
is —1 within the 21 and £ phases, and 0 within the 58 phase. 
Note, however, that 0^2 xZ 2 docs not distinguish the 21 and £ 
phases; it only implies that both are of the SPT type. Thus 
we need to look for another topological order parameter to 
distinguish these phases. 


= e 
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FIG. 5. (Color online). Topological order parameters C^ 2 xZ 2 and 
Of^ for the paths {Ti,T 2 }, (see the main text). The plot shows (i) 
both 21 and £ phases are protected by the Z 2 x Z 2 symmetry, (ii) 
the complex-conjugate symmetry only protects the 21 phase, and (hi) 
both symmetries are broken within the 23 phase. 


The ground state of the cluster Hamiltonian (19) has an ex¬ 
act iMPS form given 

f(o) = (1 - *) (^ -1 J ^ , r(i) = (1 + *) (^ J ) . 

(38) 

Since this iMPS is symmetric under the complex-conjugate 
symmetry Eq. (34) should hold. One can obtain (see 
Sec. IV A) that and = —i. Moreover, 

(cr^)((T^)* = — 1, which immediately implies — 1, 

demonstrating the SPT phase 21 is protected by Nonethe¬ 
less, we show that the £ phase is not protected by this sym¬ 
metry. 

The iMPS form of the cluster Hamiltonian (18) is expressed 
as follows: 

f (0) ^ -1 1 ^ ^ f (1) ^ 1 -1 ^ (39) 

This iMPS respects the complex-conjugate symmetry and 
Eq. (34) is obviously satisfied by = 1 and = 1. 
Hence, for this phase, = 1, implying that this phase is not 
protected by 

Summarizing, the topological order parameter is — 1, 
0, and 1 for the 21, 23, and £ phases, respectively (note that 
= 0 implies that symmetry has been broken). As de¬ 
picted in Eig. 5, has been numerically calculated through 
the paths {3i,32}- It demonstrates that the topological or¬ 
der parameter takes different values for each phase, thus 
it can truly (topologically) distinguish all three phases. This 
observation also indicates that one cannot adiabatically con¬ 
nect the £ and 21 phases because they belong to different SPT 
classes. 


V. ID COMPASS MODEL 

The ground space of the ID compass model (Eq. 2) for 
Jv = 0 can be represented as follows: 

-K = {| • • • □,□,+! • • •) I !□*) = + Al 

(40) 

where jOi) is defined on the ith blue link, and {ai,Pi} are 
two arbitrary normalization factors (|aip + |/3ip = 1). The 
ground-state subspace “Vt, with degeneracy I'Kj = 2^/2, is 
stabilized by the symmetry X (Eq. 29), that is, 

Xl'ifg) = |T-g), VlVf-g) e -K. (41) 

Thus, the symmetry group Zf = {X, 1} is respected at Jv = 
0. Turning on the coupling parameter Jv does not affect this 
fact—that the symmetry X stabilizes the ground space—up 
to the quantum phase transition point Jv = Jb (see Eig. 1- 
(d)). This behavior is due to the fact that within blue-compass 
phase the gap is not closed. However, other symmetries of the 
compass model, such as .Z = J([|^ ’ time-reversal, 

and transnational invariance, are broken here. Eor instance, Z> 
does not stabilize the ground space ‘K; rather, it rotates the 
elements of this space within itself, 

Z\ - ■ ■ niDi+i • • •) = I • • • DiDi+i • • •), (42) 

Pi) =ft|(43) 

That is, the symmetry group = {Z, 7} has been broken 
at Jv = 0. Similarly, because of the nonvanishing gap, this 
property is expected to hold within the blue-compass phase 
(see Eig. l-(d)). 

On the other hand, the symmetry Z> is stabilized by the 
ground space of the compass model at Jb = 0. This ground 
space is given by 

T7 = {|... ...) I !■) = a'l tt) + /3'l U)}, (44) 

where is defined on the ith violet link, and {«',/?'} are 
two arbitrary normalization coefficients. Eor nonzero values 
of Jb, the symmetry group is respected throughout the 
violet-compass phase as long as the gap is nonzero. In ad¬ 
dition, the symmetry group Z| rotates elements of the ground 
space it becomes broken within the violet-compass phase. 
Thus the symmetry Z 2 (Z 2 ) is preserved within the blue- 
compass (violet-compass) phase, and is broken within violet- 
compass phase (blue-compass phase). 

To verify this observation, the parameters {^Xi^z} h^^e 
also been numerically calculated using the iTEBD method, 
which are shown in Eig. 6-(a). As expected, A;^ (A^;) is equal 
to one throughout the blue (violet) phase, and is less than one 
otherwise. The breakdown of the aforementioned symmetry 
upon quantum phase transition can be captured by a local or¬ 
der parameter. Within the blue phase, where X is preserved, 
the local order parameter = (l/N) vanishes, 

(O,) = {0,X) = -{XO,) = -{O,} 

=> (Oz) = -{Oz) = 0 . 
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- 1 |) 


FIG. 6. (Color online). Numerical analysis of the symmetries 
{X, Z}, local order parameters {Oz, Ox}, and their scaling, (a) Pa¬ 
rameters {Ajf and A^;} show qualitative behavior of the symmetries 
{X and Z}, respectively. Whenever these parameters become less 
than one, the associated symmetry is broken, and when they are equal 
to one, it means the associated symmetry is preserved, (b) The local 
order parameters {Oz and Ox} vs J^/Jh- (c) Scaling of the order 
parameters near the critical point Jv/ Jh = 1. These scalings indi¬ 
cates that the magnetization exponent /3 for both order parameters is 
0.125 ± 0.001 


Whereas, the local order parameter Oz becomes nonzero 
within the violet phase, which signals the breakdown of X 
symmetry. As an example, for the ground state | • • • ft ''') 
at Jb = 0, (9^ = 1. Hence, Oz is a proper local order pa¬ 
rameter to represent the breakdown of X symmetry. In the 
same manner, one can show that the local order parameter 
Ox = ^ vanishes within the violet-compass phase 

(which comes from the preservation of the symmetry Z), and 
becomes nonzero within the blue-compass phase. We have 
numerically calculated the order parameters {Oz and Ox} by 
employing the iTEBD technique. As shown in Fig. 6-(b), 
the local order parameters {Oz, Ox} behave exactly as ex¬ 
plained before: Ox vanishes throughout the violet-compass 


phase due to the preservation of the Z symmetry, while it be¬ 
comes nonzero in the blue-compass phase as a result of Z 
symmetry breaking. Similar behavior is also observed for the 
order parameter Oz . 

Both local order parameters {Oz and Ox} vanish as the 
quantum critical point is approached JvfJh = 1, 

O, - iJv/Jb-ll'^', (45) 

Ozr^\J^/Jb-lf, (46) 

as shown in Fig. 6-(c), where /?' = 0.125 ± 0.001. It justifies 
that the corresponding order parameter exponent is /?' = 1/8, 
associated to the a second-order quantum critical point. Thus 
the local order parameters Oz and Ox faithfully capture the 
relevant physics of the ID compass model. These results are 
in agreement with Refs. 27 and 43, where physical quantities 
and phase characterization have been obtained by employing 
nonlocal string order parameters. 


VI. SUMMARY AND CONCLUSIONS 

In this paper the topological classification of the phases 
and the associated quantum phase transitions in the compass 
ladder and ID compass models have been presented by em¬ 
ploying degenerate perturbation theory, symmetry fractional- 
ization, and numerical investigation. For each phase of the 
model (denoted by 21, 25, and C), we have derived an effective 
Hamiltonian based on degenerate perturbation theory. The 21 
and £ phases have been shown to be described by two cluster 
Hamiltonians written in different bases, whereas the *8 phase 
has been shown to be represented by an Ising Hamiltonian. 
The cluster phase (specified by the ground state of the clus¬ 
ter model) is an SPT phase protected by a specific Z 2 x Z 2 
symmetry, whereby we have assigned a set of labels to specify 
them. In other words, the set of unique labels of cluster phases 
have proven to be similar to that of the 21 and £ phases. How¬ 
ever, the 21 and £ phases do not belong to the same class of 
an SPT phase: one of the phases is protected by the complex- 
conjugate symmetry, while the other is not. This observation 
has been verified by both numerical computations and analyt¬ 
ical calculations. 

We have shown that the 25 phase is of topologically-trivial 
Z 2 -symmetry breaking type, characterized by a local order 
parameter. Having determined the form of the local order 
parameter and broken symmetry, we have concluded that (i) 
the phase diagram of model is characterized by a local order 
parameter, (ii) the quantum phase transition is associated to 
a spontaneous symmetry breaking (not topological), and (iii) 
the class of the quantum phase transition is in the Ising univer¬ 
sality class, where the magnetization exponent is equal to 1 /8 
and its type is of second order. We have also verified these 
observations numerically. 

In addition, we have shown that the quantum phase tran¬ 
sition in the ID compass model (a limiting case of the com¬ 
pass ladder) is accompanied by a partial symmetry breaking. 
Fach of the phases of the ID compass model (the blue- and 
violet-compass phases) have been shown to respect part of a 
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symmetry group; the blue- and violet-compass phases respect 
two different Z 2 -symmetry groups and break the other sym¬ 
metries. Hence, upon the quantum phase transition, one of the 
Z 2 symmetries is broken and other one is preserved. This par¬ 
tial symmetry breaking has been captured by local order pa¬ 
rameters. By using numerical computations, we have shown 
that these local order parameters truly capture the quantum 
phase transition (as well as partial symmetry breaking) and 
its universality class (i.e., /3 = 1/8). It is worth mentioning 
that this type of quantum phase transition is different from, 
e.g., transverse field Ising model, in which the whole symme¬ 
try group is being broken and system goes from a disordered 
phase to an ordered one. 

Our phase characterization of the compass ladder model has 
also revealed the nature (topological classes) of a number of 


the phases of the ‘Kitaev model on arbitrary-row brick wall 
lattice’—which is similar to that of the compass ladder model. 
Although the phase diagram of this model has been known, 
the nature of remaining phases and their corresponding quan¬ 
tum phase transitions are still largely unknown. An analysis 
based on our approach, especially symmetry fractionalization 
and degenerate perturbation theory, may shed some light on 
this direction. 


ACKNOWLEDGMENTS 

This work was supported in part by Sharif University of 
Technology’s Office of Vice President for Research. 


* haghshenas@physics.sharif.edu 

' N. Goldenfeld, Lectures on Phase Transitions and the Renormal¬ 
ization Group (Perseus Books, Reading, MA, 1992). 

^ D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. Lett. 
48, 1559 (1982). 

^ C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005). 

'' N. Read and S. Sachdev, Phys. Rev. Lett. 66, 1773 (1991). 

^ A. Kitaev, Ann. Phys. 303, 2 (2003). 

® C. Nayak, S. H. Simon, A. Stem, M. Freedman, and 
S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008). 

X. G. Wen and Q. Niu, Phys. Rev. B 41, 9377 (1990). 

* X.-G. Wen, Quantum Field Theory of Many-body Systems: From 
the Origin of Sound to an Origin of Light and Electrons (Oxford 
University Press, Oxford, 2007). 

X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B 82, 155138 

( 2010 ). 

X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B 83, 035107 

( 2011 ). 

“ X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. Rev. B 84, 235128 

( 2011 ). 

N. Schuch, D. Perez-Garcia, and I. Cirac, Phys. Rev. B 84, 
165139 (2011). 

F. Pollmann, A. M. Turner, E. Berg, and M. Oshikawa, Phys. Rev. 
B 81, 064439 (2010). 

G. Vidal, Phys. Rev. Lett. 98, 070201 (2007). 

U. Schollwdck, Ann. Phys. 326, 96 (2011). 

J. Haegeman, D. Perez-Garcia, I. Cirac, and N. Schuch, Phys. 
Rev. Lett. 109, 050402 (2012). 

F. Pollmann and A. M. Turner, Phys. Rev. B 86, 125441 (2012). 
M. A. Levin and X.-G. Wen, Phys. Rev. B 71, 045110 (2005). 
C.-H. Lin and M. Levin, Phys. Rev. B 89, 195130 (2014). 

P. Fendley and E. Fradkin, Phys. Rev. B 72, 024412 (2005). 

A. Kitaev, Ann. Phys. 321, 2 (2006). 

E. K.-H. Lee, R. Schaffer, S. Bhattacharjee, and Y. B. Kim, Phys. 
Rev. B 89,045117 (2014). 

J. Chaloupka, V. Checkrelse, G. Jackeli, and G. Khaliullin, Phys. 
Rev. Lett. 110, 097204 (2013). 

J. Reuther, R. Thomale, and S. Trebst, Phys. Rev. B 84, 100406 

( 2011 ). 

J. Osorio Iregui, P. Corboz, and M. Troyer, Phys. Rev. B 90, 


195102 (2014). 

M. Barkeshli, H.-C. Jiang, R. Thomale, and X.-L. Qi, Phys. Rev. 
Lett. 114, 026401 (2015). 

X. -Y. Feng, G.-M. Zhang, and T. Xiang, Phys. Rev. Lett. 98, 
087204 (2007). 

W. Brzezicki, J. Dziarmaga, and A. M. Oles, Phys. Rev. B 75, 
134415 (2007). 

L. -M. Duan, E. Dernier, and M. D. Lukin, Phys. Rev. Lett. 91, 
090402 (2003). 

“ J. Q. You, X.-F. Shi, X. Hu, and F. Nori, Phys. Rev. B 81, 014505 

( 2010 ). 

A. Saket, S. R. Hassan, and R. Shankar, Phys. Rev. B 82, 174409 

( 2010 ). 

Y. Tserkovnyak and D. Loss, Phys. Rev. A 84, 032333 (2011). 

Y. -C. He and Y. Chen, Phys. Rev. B 88, 180402 (2013). 

F. L. Pedrocchi, S. Chesi, S. Gangadharaiah, and D. Loss, Phys. 
Rev. B 86, 205412 (2012). 

V. Karimipour, Phys. Rev. B 79, 214435 (2009). 

V. Karimipour, L. Memarzadeh, and P. Zarkeshian, Phys. Rev. A 
87, 032322 (2013). 

A. Langari, A. Mohammad-Aghaei, and R. Haghshenas, Phys. 
Rev. B 91, 024415 (2015). 

D. L. Bergman, R. Shindou, G. A. Fiete, and L. Balents, Phys. 
Rev. B 75, 094403 (2007). 

W. Son, L. Amico, R. Fazio, A. Hamma, S. Pascazio, and V. Ve- 
dral, Europhys. Lett. 95, 50001 (2011). 

D. V. Else, I. Schwarz, S. D. Bartlett, and A. C. Doherty, Phys. 
Rev. Lett. 108, 240505 (2012). 

S. Montes and A. Hamma, Phys. Rev. E 86, 021101 (2012). 

Z. Nussinov and J. van den Brink, Rev. Mod. Phys. 87, 1 (2015). 

H. T. Wang and S. Y. Cho, J. Phys.: Condens. Matter 27, 015603 
(2015). 

M. Kargarian, H. Bombin, and M. Martin-Delgado, New J. Phys. 
12, 025018 (2010). 

M. B. Hastings, J. Stat. Mech.: Theo. Exp. , P08024 (2007). 

R. Orus and G. Vidal, Phys. Rev. B 78, 155117 (2008). 

R. Haghshenas, A. Langari, and A. T. Rezakhani, J. Phys.: Con¬ 
dens. Matter 26, 456001 (2014). 

D. Perez-Garcia, F. Verstraete, M. M. Wolf, and J. I. Cirac, Quan¬ 
tum Inf. Comput. 7, 401 (2007). 



